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We propose an adjoint SU (5) GUT model with a Ty family symmetry and an extra ^2 ® ^3 < 8 ) -^4 < 8 ) 

•^4 < 8)^12 discrete group, that successfully describes the prevailing Standard Model fermion mass and 
mixing pattern. The observed hierarchy of the charged fermion masses and the quark mixing angles 
arises from the Za® . 248)^12 symmetry breaking, which occurs near the GUT scale. The light active 
neutrino masses are generated by type-I and type-III seesaw mechanisms mediated by the fermionic 
SU{5) singlet and the adjoint 24-plet. The model predicts the effective Majorana neutrino mass 
parameter of neutrinoless double beta decay to be mpp = 4 and 50 meV for the normal and the 
inverted neutrino spectra, respectively. We construct several benchmark scenarios, which lead to 
SU{5) gauge coupling unihcation and are compatible with the known phenomenological constraints 
originating from the lightness of neutrinos, proton decay, dark matter etc. These scenarios contain 
TEV-scale colored helds, which could give rise to a visible signal or be stringently constrained at 
the LHG. 


I. INTRODUCTION 


The LHC discovery of a 126 GeV Higgs boson [I] has completed the era of the experimental quest for the missing 
elements of the Standard Model (SM). Now it is a well- established and extremely successful theory of the electroweak 
phenomena. However, the SM has several unaddressed issues, such as the smallness of neutrino masses, the observed 
pattern of fermion masses and mixings, and the existence of three fermion families [2], etc. The observed pattern of 
fermion masses spans over a range of 5 orders of magnitude in the quark sector and much wider considering neutrinos. 
Neutrino oscillation experiments demonstrated that at least two of the neutrinos are massive, with masses much lower, 
by several orders of magnitude, than the other SM fermions, and also that all of the three neutrino flavors mix with 
each other. The smallness of quark mixing contrasts with the sizable mixing of neutrinos; i.e., while in the quark 
sector all the mixing angles are small, in the neutrino sector two of them are large, and only one mixing angle is small. 
This suggests that the neutrino sector is described by a different kind of underlying physics than the quark sector. As 
is well known the tiny neutrino masses may point towards a high-energy scale of New Physics, where lepton number 
violation (LNV) takes place. 

The physical observables in the neutrino sector, i.e., the neutrino mass-squared splittings and mixing parameters, 
are constrained from the global fits of the available data from neutrino oscillation experiments Daya Bay [S], T2K 
[4], MINOS [5], Double CHOOZ [6], and RENO [7], as shown in Tables |l] and | h] (based on Ref. [8]) for the normal 
(NH) and inverted (IH) hierarchies of the neutrino mass spectrum. As seen from these tables the neutrino oscillation 
experimental data show a clear evidence of a violation of the so-called tribimaximal symmetry described by the 
tribimaximal mixing matrix (TBM), predicting the neutrino mixing angles (sin^ 012) (sin^023)y^^ = 5, 

and (sin^ ^i3)tbm = O' 
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Parameter 

Amii( 10 -®eV 2 ) 

Amii(10-3eV2) 

(sin^ 6 * 12 ) 

V / exp 

(sin^ 623) 

V /exp 

(sin^ 013 ) 

V /exp 

Best fit 

7.60 

2.48 

0.323 

0.567 

0.0234 

Icr range 

7.42 - 7.79 

2.41 - 2.53 

0.307 - 0.339 

0.439 - 0.599 

0.0214- 0.0254 

2 a range 

7.26 - 7.99 

2.35 - 2.59 

0.292 - 0.357 

0.413-0.623 

0.0195 - 0.0274 

3a range 

7.11 - 8.11 

2.30 - 2.65 

0.278 - 0.375 

0.392 - 0.643 

0.0183 - 0.0297 


Table I: Range for experimental values of neutrino mass squared splittings and leptonic mixing parameters, taken from Ref. 
[S], for the case of normal hierarchy. 


Parameter 

Amii( 10 -®eV^) 

Am? 3 ( 10 -®eV 2 ) 

(sin^ 612 ) 

V / exp 

(sin2 023)^,„ 

^13)exo 

Best fit 

7.60 

2.38 

0.323 

0.573 

0.0240 

la range 

7.42 - 7.79 

2.32 - 2.43 

0.307 - 0.339 

0.530 - 0.598 

0.0221 - 0.0259 

2 a range 

7.26 - 7.99 

2.26 - 2.48 

0.292 - 0.357 

0.432 - 0.621 

0.0202 - 0.0278 

3a range 

7.11 - 8.11 

2.20 - 2.54 

0.278 - 0.375 

0.403 - 0.640 

0.0183 - 0.0297 


Table II: Range for experimental values of neutrino mass squared splittings and leptonic mixing parameters, taken from Ref. 
[S], for the case of inverted hierarchy. 

Addressing the flavor puzzle requires extensions of the SM, including larger scalar and/or fermion sectors, as well as 
an extended gauge group with additional flavor symmetries, which allow one to explain the SM fermion mass and 
mixing pattern. Along this line, several models have been proposed in the literature (for a review see, e.g.. Refs, 
[nn]). The fermion mass and mixing pattern can also be described by postulating particular mass matrix textures 
(see, e.g., Ref. [14] for a comprehensive review and some recent works considering textures). It is believed that grand 
unified theories (GUTs) endowed with flavor symmetries could provide an unified description for the mass and mixing 
pattern of leptons and quarks. This is motivated by the fact that leptons and quarks belong to the same multiplets 
of the GUT group, allowing to relate their masses and mixings [HHII]- Furthermore, this setup can generate small 
neutrino masses through a type-I seesaw mechanism, where the new heavy Majorana neutrinos acquire very large 
masses due to their interactions with scalar singlets, assumed to get vacuum expectation values (VEVs) at the very 
high energy scale. Various GUT models with flavor symmetries have been proposed in the literature [I^37j . For a 
general review, see for example [351 EH] ■ 

In this paper we propose a version of the adjoint SU{5) GUT model, where an extra ^2 0 Z 3 0 Z 4 ® Z 4 0 Z 12 
discrete group (see Ref [40] for a comprehensive study of the Ty flavor group) extends the symmetry of the model, and 
several scalar fields are included to generate viable and predictive textures for the fermion sector. The particular role 
of each additional scalar field and the corresponding particle assignments under the symmetry group of the model are 
explained in detail in Sec. |llj The fermionic sector of our model, in addition to the SM fermions, contains one heavy 
Majorana neutrino Nn, singlet under the SM group, and an adjoint 24 fermionic irreducible representation (irrep) of 
SU{5) so that the light neutrino masses are generated via type-I and type-III seesaw mechanisms. An adjoint SU{5) 
GUT model, but without discrete symmetries, and more minimal particle content was previously considered in Ref. 
[4T] . Although our model is less minimal than that of Ref. [41] , it provides a successful description of the SM fermion 
mass and mixing pattern, not addressed in Ref. |4Ij. Also, our model is more predictive in the leptonic sector. 

The model has 14 free effective parameters, which allow us to reproduce the experimental values of 18 observables with 
good accuracy: i.e., 9 charged fermion masses, 2 neutrino mass-squared splittings, 3 lepton mixing parameters, and 4 
parameters of the Wolfenstein CKM quark mixing matrix parametrization. It is noteworthy that the alternative SU (5) 
GUT model of Ref. [25] . with a supersymmetric setup and flavor symmetries, also has 14 free effective parameters 
aimed at reproducing the above mentioned 18 observables. Discrete symmetries different from Ty have also been 
employed in adjoint SU{5) GUT models. Some examples are the supersymmetric adjoint SU{5) GUT model with 
A 4 flavor symmetry of Ref. [33], and the nonsupersymmetric SU{5) GUT model with Z 4 symmetry of Ref. [30j . 
The aforementioned SUSY adjoint SU{5) GUT model with A 4 flavor symmetry employs two sets of Z 2 symmetries 
and two sets of U{1) symmetries. One of these U{1) symmetries is global and represents an i?-parity symmetry, 
whereas the other C/(l) symmetry is assumed to be gauged. The two Z 2 symmetries shape the Yukawa matrices 
for quarks and charged leptons. That model includes Higgs multiplets in 1 , 5 , 5 , 45 , dimensional representations 
of SU{5). The 14 SU{5) scalar singlets are grouped into three A 4 triplets and two A 4 trivial singlets. As in our 
model, neutrino masses arise from a combination of type-I and type-III seesaw mechanisms. Whereas in the model 
of Ref. |29j . the CKM quark mixing matrix mainly arises from the down-type quark sector, in our model, the quark 
mixing is completely determined from the up-type quark sector. Furthermore, in this model the leptonic mixing 
angles are determined by two parameters, whereas in our model only one parameter determines the leptonic mixing 
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angles. Moreover, the quark masses and mixings are studied in detail in our model, whereas in the model of Ref. [55] 
a detailed study of quark masses and mixings is not included. Besides that, our model includes a detailed discussion 
of gauge coupling unification and seesaw mass scale limits, not performed in the model of Ref. [29]. With respect to 
the nonsupersymmetric adjoint SU{5) GUT model with Z 4 discrete symmetry of Ref. |30j . the discrete symmetry is 
introduced in that model in order to generate the nearest-neighbour-interaction textures for charged fermions. The 
scalar sector of that model includes an adjoint multiplet, a quintuplet and one 45-dimensional representation, and 
neutrino masses arise from type-I seesaw, type-III seesaw, and one-loop radiative seesaw mechanisms. Despite the 
fact that this model is more minimal than our T 7 flavor adjoint SU{5) GUT model, our model has a much more 
predictive lepton sector. Let us note that the lepton sector of the model of Ref. [30] has 12 effective free parameters, 
whereas the lepton sector of our model has a total of 6 effective free parameters. Unlike the models of Refs. [29l I30j . 
our adjoint SU{5) GUT model, which is non supersymmetric, employs a Z 2 G Z 3 0 ^4 0 Z 4 0 Z 12 discrete symmetry. 

After imposing gauge coupling unification equal to or better than that in the MSSM (another difference with respect 
to Refs. [311142], which impose exact gauge coupling unification), we find a wide set of simple field configurations, 
which pass proton decay constraints and give rise to the neutrino masses through a type-I and type-III seesaw 
realization. Considering the limit on the triplet scalar mass (denoted here as S 3 ), which comes from the cold dark 
matter constraints pointed out in Ref. [43] . we found lower limits on the seesaw scale for two simple scenarios with 
different sets of beyond-the-SM field configurations. 

The paper is organized as follows. In Sec. [Tl]we outline the proposed model. In Sec. |III| we present our results on 
neutrino masses and mixing, followed by a numerical analysis. Our results for the quark sector, with the corresponding 
numerical analysis, are presented in Sec. [IV| Gauge coupling unification and seesaw scale mass limits are discussed in 
Sec[V] We conclude with discussions and a summary in Sec. m Some necessary facts about the Tj group and details 
of our analysis are collected in appendixes. 


II. THE MODEL 


The first grand unified theory (GUT) proposed in Ref. [33] is based on the SU{5) gauge symmetry accommodating 
the SM fermions in 5 -I- 10 and the scalars in 5 -|- 24 irreps of SU{5). As is well known this model suffers from 
several problems. In particular, it predicts wrong down-type quark and charged lepton mass relations and a short 
proton lifetime, and the unification of gauge couplings disagrees with the values of as, sin^vv and aem measured at 
the Mz scale. Moreover this minimal SU{5) GUT model does not include a mechanism for generating nonvanishing 
neutrino masses in contradiction with experimental data on neutrino oscillations. The minimal iS't7(5) GUT model 
can be improved by including, in particular, a scalar 45 irrep of SU{5) [42l l43l I45H55] . This next-to-minimal SU{5) 
GUT model fails, however, in describing the observed pattern of fermion masses and mixings, due to the lack of 
explanation for the hierarchy among the large number of Yukawa couplings in the model. Below we consider a multi- 
Higgs extension of the next-to-minimal SU (5) GUT model, which successfully describes the pattern of the SM fermion 
masses and mixing. The full symmetry Q of the model is broken in two subsequent steps as follows: 

g = SU (5) (g) Ty (g) Z 2 G Z 3 (g) Z 4 (g) (g) Z 12 (I) 

'll ^GUT 

SU {3)cZ>SU{2)^Z>U{1)y 
'll ^EW 


^C/(3)cGC/(I)_ 


Let us note that among the discrete symmetries we introduced the non-Abelian flavor symmetry group Ty, which is 
the smallest group with a complex triplet representation, allowing us to naturally accommodate the three families of 
fermions. 

The fermion assignments under the group G = SU{5) ® Z^® Z^® Z^® Z'^® Z 12 are 


4- 


( 1 ) 


(10, lo,I,a;, I, l,i) 


4 ^ 


( 2 ) 


( 


10,ll,I,a;^I,I,e- 


4 - 


(3) 


(10,12,1,1,1,1,1), i,j = l,2,3,4,5 (2) 


^ , iVyy-(l,lo,-1,1,1,-1,-1), p^(24,lo,-l,l,l,-i,l), 


where lo = 


2TTi 

e 3 . 
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More explicitly, the fermions are accommodated as 



/ 

0 





\ 

1 

71 



0 



- 4 ') 





0 


-4^) 




4^^ 

4 ^) 

0 

-lifU 



1 


4^) 

4 ^^ 

lifu 

0 

/ 


/=1,2,3 z,j = l,2,3,4,5. (3) 




(4) 


Here the subscripts correspond to the different quark colors, while the superscript / refers to fermion families. One 
can see that the three families of left- and right-handed fermions, corresponding to the 5 irrep of SU{5), are unified 
into a T 7 antitriplet 3, while the three families of left- and right-handed fermions corresponding to the 10 irreps of 
SU (5) are assigned to the three different Tj singlets lo, li, l 2 - 

The scalar sector is composed of the following SU (5) representations: one 24, one 45, four 5’s and ten I’s. Two sets 
of SU (5) singlets are unified into two Ty triplets. The remaining scalar fields, i.e., one 45, one 24, four 5’s and the 
remaining four I’s, are accommodated by three Ty singlets. Thus the Q assignments of the scalar fields are 


a 


hW 

^(4) 


(l,lo,l,l,l,l,e-T) , 

(l,lo,l,a;,- 1 , 1 , 1 ), 

(5,lo,-l,l,l,l,e-^) 

(5,12,1,1,1,1,1), 


(l,lo,l,a;,f,l,e-Tj , ^ ^ 2 , ^ 3 )(1> 3,-1,1,1,1,1), 

(1,1,1,1,1,-1,1), X = (Xi,X 2 ,X 3 ) ^ (l,3,l,l,l,f,e'?) , 
ilf) ^ (5, lo,l, w, 1,1,1), fP ~ (5, li,l,1,1, l) , 

E] ~ (24, lo,l, 1,1,1,1) , ^ ( 45 , lo, -1,1,1,1, e-^) . 


(5) 


The VEVs of the scalars (h = 1, 2, 3,4) and S® are 



(Ux 


h = l,2,3,4. 






z,j = l,2,3,4,5. (6) 


Note that the VEV pattern for the S held given above is consistent with the minimization conditions of the model 
scalar potential and follows from the general group theory analysis of spontaneous symmetry breakdown |56j . 

The following comments about the possible VEV patterns for the Ty scalar triplets x and ^ are in order. Here we 
assume a hierarchy between the VEVs of the Ty scalar triplets x and i.e., << , which implies that the mixing 

angle between the the Ty scalar triplets x and ^ is strongly suppressed since it is of the order of (N, as follows from 
the method of recursive expansion of Refs. |57H59] . Consequently, the mixing between the Ty scalar triplets x and ^ 
can be neglected. The parts of the scalar potential for each of the two Ty scalar triplets at the renormalizable level 
are given by: 


rW ^ _^2 (xx*)3 (xx*)3 + %.2 (xx)3 (x*x*)3 + %.3 (XX)3 (X*X*)3 

+KxA (xx*)i0 (xx*)i0 + %.5 (xx*)ii (xx*)i2 + H.c. ( 7 ) 

^7^ = -4 (er)i„ + (cr )3 in+ms icth+icOs m^ 

+««.4 (?r)i„ m)io + ««.5 (cDi, (cDi, + H.c. (8) 


In the part of the scalar potential for each scalar triplet there are six free parameters: one bilinear and five quartic 
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couplings. The minimization conditions of and 


T7 


lead to the following relations: 


dvsi 


dvs. 


dvss 


- 2 nsi ^5 + {vs 2 + uh) + 4 (ks ,2 + ks.s) cos (6<Si - 0s^) + Vg^ cos (6»Si - 6»sJ] 

+8ks,4VSi {vg^ + Vg^ + Vg^) +4ks,5'CSi {‘^Vg^-Vg^ - VgJ 
0 

- 2 vg^Hg + ‘inSAVg^ {vg^ + VgJ + 4 {Kg ^2 + KS.s) l^Si COS (05, - 05J + Vg^ COS (052 - 053)] 
+8KSAVS2 {vs, + Vg, + UI3) +4k5,5^^S2 (242- 4 i - ^k) 

0, (9) 

- 2 vs^fis + 4 K 5 ,in 53 (u|, + n|J + 4 (55,2 + ks.s) [ 3 i^S 3 'cs 2 cos (05^ - 9 g^) + cos (05, - 053)] 

+8k5,4^^S3 {vk + + ^k) +4«S,5?^S3 ( 243 -4i - 

0 , 


where m = 1,2, S = XjC and (S') = . Then, from an analysis of the minimization 

equations given by Eq. ([^ and setting obtain for a large range of the parameter space the following 

VEV direction: 


^Xl 


= p. 2 


t!. 


X 

72’ 


V 


^X3 


X 

72’ 


‘"X2 


= 0 , 


= ^2 = = 


73 - 


( 10 ) 


In the case of this is a vacuum configuration preserving a Z 3 subgroup of T-j, which has been extensively studied 
by many authors (see, for example. Ref. |40j 1. The VEV pattern for x is similar to the one we previously studied 
in an SU{5) model and in a 6HDM with A 4 flavor symmetry [331 US]- It is worth mentioning that there could be 
relative phases between the different components of (^), consistent with the scalar potential minimization equations, 
as follows from the expressions given in Eq. We have checked that the nonvanishing phases consistent with the 
scalar potential minimization equations satisfy Og. = 7 ^Ski with i ^ j ^ k (i,j,k = 1,2,3). Moreover, we 

checked that the physical observables in the lepton sector, studied in Sec. El do not depend on these phases. For 
generality we included nonzero phases in the (y) sector, as indicated in Eq. (10). 

From the expressions given in Eq. (§, and using the vacuum configuration for the Tj scalar triplets given in Eq. ( jlOj ) , 
we find the relation between the parameters and the magnitude of the VEV: 


9^1 = 


«x.i 


■ 4k 


xA ' 


''X.Sj 


+ 2 (k5,2 + ^^{, 3 ) + 3kj,4] 


( 11 ) 


These results show that the VEV directions for the Ty triplets y and ^ in Eq. 
minimum of the scalar potential of our model. 


(10) are consistent with a global 


Assuming that the charged fermion mass pattern and quark mixing hierarchy is caused by the Z 3 , Z 4 , and Z 12 
symmetries, and in order to relate the quark masses with the quark mixing parameters, we set the VEVs of the SU{5) 
scalar singlets as follows: 


= Vr = = Va = AquT = AA, 


( 12 ) 


where A = 0.225 is one of the parameters in the Wolfenstein parametrization and A is the high-energy scale cutoff 
of our model, to be clarified below. Assuming that the parameters of the scalar interaction terms involving these 
SU{5) scalar singlets are of the same order of magnitude, it is straightforward to show that the VEVs in Eq. (12) are 
consistent with the minimization conditions of the model scalar potential. 

The fields being the 45 irrep of SU{5), satisfy the following relations: 


®;i. = -*L,. 


E*t = o. 


2 = 1 


Consequently, the only allowed nonzero VEVs of are 


hj,k = 1,2, - • • ,5. 


(13) 


(‘I’ps) — q 


{S)-AS\S^) 


i,j = 1,2, 3,4, 5, 


p = l,2,3,5. (14) 
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With the specified particle content, there are the following interaction terms, invariant under the group 1 
for the further analysis: 

r o /i o A 


Ly 




+ft (</■•{).. ^ 

r /1\ 3 , 313 (T® 313 3 3 303 3T® 303 3 , 31 3 iT^ 


I jk ^3 


A; 




^ ^23 

[V'* (X*X*)3] 1 „ ^ 

A413 




L- V.X- .V ^OJlQ 6 l-\- /J-- J iQ «- yY ^ '■' i-Q J 

(V'*x)i„ ‘^%P^k + rriNNuNfi + yiN rN'^ '^ a + xit t + X2ip >p ^ (15) 


where the dimensionless couplings in Eq. (15) are 0(1) parameters, and are assumed to be real, excepting 7 ^ 2 , 
721 , 731 and 733 , which are assumed to be complex. The subscripts Iq, li, I 2 denote projecting out the corresponding 
T-j singlet in the product of the two triplets. Let us note that Eq. (15) is SU (5) invariant, since the scalar fields 
U(h)p _ 1 ^ 2 , 3,4), $*^., transform as 5, 5, 45 and 45 under SU (5), respectively, and the fermionic fields '0* 

f _ 

and (/ = 1,2,3) as 5 and 10 , under the SU (5) group, respectively. Besides that, it is worth mentioning that 

the scalar field (h — 1, 2, 3,4) transforms with the opposite Zri charges as compared to Furthermore, it is 

noteworthy that the term Tr (p^) is not present in Eq. (15), since this term is not invariant under the Z'^ symmetry. 
The lightest of the physical neutral scalar states of and 4>, should be interpreted as the SM-like 

126 GeV Higgs observed at the LHC. 


Let us summarize and comment on the above presented model setup. In comparison with the next-to-minimal SU{5) 
GUT model of Refs. besides for the introduction of additional discrete symmetries, we also extended 

the fermionic sector by introducing one heavy Majorana neutrino Nr, singlet under the SM group and a 24 fermionic 
irrep of SU (5), namely pj. We will show that since the Z 2 symmetry present in (j^ is not preserved at low energies, 
the active neutrinos get tree-level masses via type-I and type-HI seesaw mechanisms. In the next section, we will 
also show that in order to successfully accommodate the experimental data on neutrino mass-squared splittings, one 
needs both the SM singlet right-handed Majorana neutrino and the 24 fermionic irrep of SU{5). Having only one of 
them would lead to two massless active neutrinos, in contradiction to the experimental data on neutrino oscillations. 
Note that our fermionic sector is less minimal than the one considered in Ref. m with only a 24 fermionic irrep of 
SU (5). However, our model provides a successful description of the SM charged fermion masses and mixing pattern, 
not addressed in Ref. mi- 

Despite the flavor-discrete groups in Eq. ([^, the corresponding field assignment as well as the VEV pattern look 
rather sophisticated, although each introduced element plays its own role in the arrangement of the desired particle 
spectrum and flavor mixing. Let as briefly sketch our justification of the model setup: 


1. The scalar sector includes the following SU (5) representations: one 24, one 45, four 5’s and ten I’s. The 45 
and the four 5’s scalar irreps of SU (5) acquire VEVs at the electroweak scale, thus inducing the second step 
of symmetry breaking. The remaining scalars acquire VEVs at the GUT scale and trigger the first step of 
symmetry breaking. As previously mentioned, having scalar fields in the 45 representation of SU (5) is crucial 
to get the correct mass relations of down-type quarks and charged leptons. 

2. The T^ discrete group is crucial to generating textures for the lepton sector that successfully account for the 
experimentally observed deviation from the trimaximal mixing pattern that attracted a lot of attention in the 
literature as a framework for describing the lepton mixings; see for example Ref. |40) . To reproduce the nontrivial 
quark mixing consistent with experimental data, the up-type quark sector requires three 5’s, i.e., iL) 7 iL) 
and irreps of 51/(5) assigned to different Tj singlets. In the down-type quark sector, on the other hand, 
only one 5 irrep one 45 irrep 4)*^ assigned to T^ trivial singlets and three I’s, unified in the T 7 triplet 
are needed. 
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3. The Zi symmetry separates the scalars in the 5 and 45 irreps of SU (5) participating in the Yukawa interactions 
for charged leptons and down-type quarks from those ones participating in the Yukawa interactions for up-type 
quarks. This implies that the SU (5) scalar multiplets contributing to the masses of the down-type quarks 
and charged leptons are different from those that provide masses to the up-type quarks. Furthermore, the 
symmetry separates the scalar triplet ^ participating in the Yukawa interactions for charged leptons and 
down-type quarks, from that one (y) participating in the neutrino Yukawa interactions. In the scalar sector, 
the Z 2 symmetry distinguishes the Ty scalar triplet the SU{5) multiplets and charged under this 
symmetry, from the remaining scalar fields, neutral under this symmetry. Because of this, the 45 and one of 
the 5’s scalars participate in the Yukawa interactions for leptons and down-type quarks, whereas the remaining 
SU (5) multiplets participate in the Yukawa interactions for up-type quarks. This results in a reduction of 
parameters in the quark sector, since due to the Z 2 symmetry the 45 scalar irrep of SU (5) does not appear in 
the up-type quark Yukawa terms. Furthermore, all fermions are Z 2 even excepting the right-handed Majorana 
neutrino and the 24 fermionic irrep of SU (5), which are Z 2 odd. 

4. As with the Ty symmetry, the Z'^ symmetry is also necessary to get a predictive neutrino mass-matrix texture 
that only depends on three effective parameters and that gives rise to the experimentally observed deviation 
from the trimaximal mixing pattern. This symmetry also separates the 24 fermionic irrep p and the charged 
under this symmetry from the remaining fermionic fields, neutral under this symmetry. 

5. The Z^ and Z 4 symmetries are crucial to getting the right pattern of charged lepton and down-type quark masses. 
The Z 3 symmetry distinguishes the three lO’s irreps of SU{5), having different Z 3 charges. The Z 4 symmetry 
separates the SU (5) scalar singlets ip and r, charged under this symmetry, from the remaining scalar helds, 
neutral under this symmetry. All the fermionic fields are assumed to transform trivially under the Z 4 symmetry. 
Without the Z^ and Z 4 symmetries, the down quark and electron masses would be larger by about 2 orders of 
magnitude than their corresponding experimental values, unless one sets the corresponding Yukawa couplings 
unnaturally small. It is noteworthy that, unlike in the up-type quark sector, a A® suppression (A = 0.225 is one 
of the Wolfenstein parameters) in the 11 entry of the mass matrices for down-type quarks and charged leptons 
is required to naturally explain the smallness of the down quark and electron masses. The Z 3 , Z 4 , and Z 12 
symmetries will be crucial to achieve that A® suppression. 

6 . The Z 12 symmetry shapes the hierarchical structure of the quark mass matrices necessary to get a realistic 
pattern of quark masses and mixings. Besides that, the charged lepton mass hierarchy also arises from the Z 12 
symmetry. Let us recall that due to the properties of the Zn groups, it follows that Z 12 is the lowest cyclic 
symmetry that allows building the dimension-ten up-type quark Yukawa term with a (T®/A® insertion in a term 
of dimension four, crucial to getting the required A® suppression in the 11 entry of the up-type quark mass 
matrix. 


III. LEPTON MASSES AND MIXING 


The charged lepton mass matrix is derived from Eq. (151 by using the product rules for the Ty group given in 
Appendix]^ and considering that the VEV pattern of the SU{5) singlet T-j scalar triplet ^ satisfies Eq. (10) with the 
VEVs of their components set to be equal to AA (A being the cutoff of our model), as indicated by Eq. (12). Then, 
the mass matrix for charged leptons takes the form 


M, - 


0 

0 


0 

,( 0 x 5 



= ^iL^iag (me 




^r) 


(16) 


with 


(0 

a = ■ 


ViL = -7^ 



(0 ^ 
aU = - 


( (1) 

\a2VH 


6/32?^<e>) 


(0 

a\ = - 


^ - 6/?3U<[,) , 


(17) 


(18) 


Here A = 0.225 is the Wolfenstein parameter. As was commented in the previous section, we assume that the 
dimensionless couplings ai and /3j (* = 1,2, 3) in Eq. (15) are roughly of the same order of magnitude and the VEVs 









and Vis, are of the order of the electroweak scale v ~ 246 GeV. Therefore, the hierarchy among the charged lepton 


masses arises from the breaking of the Z 3 , Z 4 and Z 12 symmetries. As seen, the lepton mass matrix Eq. (16) is 
fully determined in our model by three effective parameters ai \ 3 shown in Eq. |l7[ ), which we fit to reproduce the 
experimentally measured values of lepton masses and mixings at the Mz scale. A similar situation takes place in the 
sector of quarks, as will be shown in the next section. 


From the neutrino Yukawa terms of Eq. (15), and taking into account that the VEVs of the SU{5) singlets ip, a and 


T are set to be equal to AA (where A is our model cutoff), as indicated by Eq. (12), we find that the fields contained 
in the 24 fermionic irrep of SU{5) acquire very large masses, which are given by 


"Po 


"P 3 


"Ps 




( 3 , 2 ) 


y2v-p - 

ysv-Vn 

yMA ’ 


3y3V~v^ 

y2v-p - 

v^A ’ 



y2v-p + 

V^A ’ 

^( 3 , 2 ) 

= y2Vr, - 


ysVsV'n 


(19) 


Here irin^, mn„ and 


ipg are the masses of the fermionic singlet Pq 

and 


triplet P 3 and octet pg contained in the 24 
the masses of the (3, 2) and (3, 2) fermionic 


"Po’ "Ps 

fermionic irrep of SU (5), respectively. We denote by 
fields corresponding to the SU{3) triplet and SU{3) antitriplet, SU{2) doublet parts of p, respectively. Consequently, 
the light active neutrino masses arise from type-I and type-III seesaw mechanisms induced by the SU{5) singlet heavy 
Majorana neutrino Nn, the fermionic singlet pg and the fermionic triplet P 3 , respectively. 


Eq. (10), we find the neutrino mass matrix: 


From the neutrino Yukawa terms of Eq. (15) and the VEV pattern of the SU{5) singlet T 7 scalar triplet x given by 


= 


O: 

(M, 


3x3 

D\'^ 



( ° 

i4> 

Eie 2 

id) 

E2e 2 

id) id) 

Y 2 e~^ Y 2 e~^ 

=\ X 

0 

0 

0 0 

V 0 

i<t> 

EieT 

i <b 

¥ 26^3 

id) id) 

¥ 26 ^ ¥ 26 ^ 


Mu = 


( rriN 
0 
0 
0 

V 0 


0 

m 


Po 

0 m 


'P 3 


0 

0 

0 

«P3 


0 \ 
0 
0 
0 

^P3 


( 20 ) 


where: 


A = (Ai, + A2.)uW^, 


El = 


3 /T 5 




A4i/U<j> 


A’ 


Eo — ( A: 


- 3A4,.U4,^ (21) 


Therefore, the light neutrino mass matrix takes the following form: 

/ 0 A 

OHO 
y A 0 Ae*'^ 




( 22 ) 


where 


3 SEl 

A = 


B = 


mjv' 


(23) 


The smallness of neutrino masses in our model is the consequence of their inverse scaling with respect to the large 
masses of the singlet pg and the triplet pg fermionic helds and proportionality to the squared neutrino Yukawa 
couplings. 
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The mass matrix in Eq. (22 1 for light active neutrinos is diagonalized by a unitary rotation matrix Vi,. There are 


two solutions of this diagonalization problem: 


/ mi 

0 

® 1 

1 1 

( cos 6 

0 

sin0e \ 

0 

m 2 

0 

, with 14 = 

0 

1 

0 

V 0 

0 

m 3 ) 

1 1 

1 -sin0e*'^ 

0 

COS0 / 


Pu 


-4 


Pu = 


diag ^ 


gi«l/2 gia2/2 ^ 1013/2 


)■ 


(24) 

(25) 


The solutions corresponding to 6 * = + 7 r /4 and 6 = — 7 r /4 we identify with the normal (NH) and inverted (IH) neutrino 
mass hierarchies, respectively, so that 


NH : e = +- : 

4 


IH : e = -: 

4 


= 0 , 

m^^ = B, 

TTiiy^ — 

ai = Q!2 

= 0 , 

as = </>, 

(26) 

= 2 H, 

mv2 = B, 

rrijy^ = 0 , 

Q!2 = 0:3 

= 0 , 

«! = -4>- 

(27) 


Let us note the presence of nonvanishing Majorana phases (j) and —cj) for NH and IH cases, respectively. This simple 
relation requires the effective dimensionful parameters A and B, given by Eq. (23), to be real, which is consistent 
with our previously mentioned assumption that the dimensionless couplings in Eq. (15) are 0(1) parameters assumed 
to be real, excepting jn, 721 ) Tsi 75 ^ 3 , assumed to be complex 
Now we find the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) leptonic mixing matrix. 


/ 


[/ = ^IK = 


cos 0 

73 


V3 


1 

V3 


COS ^ 


cos 9 
\ x/3 


3 sin 0 e 3 






cos 9 
V3 


V3 


le 




+ 


2z7r 2z7r _ 2t7r 

3 sin 0 e 3 e 3 cos i 


^/3 


V3 


^/3 


V3 

sin 9 


e sin 6 

V3 / 


4 


(28) 


with the pattern of the trimaximal (TM 2 ) type [SI]. It is noteworthy that in our model the PMNS matrix depends 
on a single parameter 0, and the neutrino masses (26) and (27) depend on two parameters, A and B. Comparing 


the matrix U in Eq. (28) with the standard parameterization of the PMNS matrix in terms of the solar 612 , the 


atmospheric 623 and the reactor 613 angles, we find 


Si4 012 = 


142 r 


l-\Ue 3 \ 2- Z -cose 


si4 013 = |t/e3|^ = 1(1 + Z • COS0), 


sin 023 = 


i43r 


2 — z • (cos (j) + 's/S sin (j)) 


l-|C/e3r 


4 — 2z • cos ( 


(29) 


with z = 1 and z = —1 for NH and IH, respectively. Note that in the limit (j> = 0 and (/) = tt for IH and NH, 
respectively, the mixing matrix in Eq. (28) reduces to the tribimaximal mixing pattern, which yields a vanishing 
reactor mixing angle 0 i 3 . 

For the Jarlskog invariant and the CP-violating phase [5], we find 


J = Im ( 41424 * 24 * 1 ) = ^ 


sin (5 = 


8J 


cos 013 sin 2012 sin 2023 sin 20 i 3 


= 0 


(30) 


since 0 = ±? according to Eq. (24). 


Thus, our model predicts a vanishing leptonic Dirac CP-violating phase. 

In what follows, we adjust the three free effective parameters (/>, A and B of the active neutrino sector of our model 
to reproduce the experimental values of three leptonic mixing parameters and two neutrino mass-squared splittings, 
reported in Tables jlj and jlTl for the normal and inverted neutrino mass hierarchies, respectively. We fit the parameter 
(j) to adjust the experimental values of the leptonic mixing parameters sin^ Oij, whereas A and B are fixed so that the 
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measured mass-s qua red differences are reproduced for the normal (NH) and inverted (IH) neutrino mass hierarchies. 
From Eqs. (27), (26), and the definition = mj — m'j, we find 


NH : iJivi = 0, TOjyj = B = ^ 

« 9meV, 

= 2 A = 1 

J Aml^ ? 

» 50meV; 

(31) 

IH : 171^2 = B = ■\J -I- « 

50meV, 

m,2i = 2 |A| = - 

\J Ato?3 ' 

« 49meV, = 0, 

(32) 

for the best-fit values of Arnfj taken from Tables [ 
To fit the leptonic mixing parameters sin^ 0^ 

i|and[H| 

to their experimental values, given 

in Tables [Ij [nj we vary 

the (j) 

parameter, finding the following best fit result: 

NH : (j) = —0.887r, sin^ 0i2 

« 0.34, 

sin^ 023 « 0.61, 

sin^ 013 

« 0.0232; 

(33) 

IH : (j)= 0.12 TT, sin2 0 

12 « 0.34, 

sin^ 023 « 0.61, 

sin^ 0 

'13 « 0.0238. 

(34) 


Thus, sin^ 9 i3 is in excellent agreement with the experimental data, for both normal and inverted neutrino mass 
hierarchies, whereas sin^ 9 i2 and sin^ 023 exhibit a 2cr deviation from their best-fit values. 

Now we are ready to make a prediction for the neutrinoless double beta {0v(3/3) decay amplitude, which is proportional 
to the effective Majorana neutrino mass parameter. 


mi3i3 


k 


(35) 


where and are the PMNS mixing matrix elements and the Majorana neutrino masses, respectively. Using 
Eqs. (24)-(28) and (31)-(34), we get for both normal and inverted hierarchies 


1 

mpp = - 


+ 4:A cos^ 



4 meV for NH 
50 meV for IH 


(36) 


These values are beyond the reach of the present and forthcoming QufijS decay experiments. The presently best upper 
limit on the effective neutrino mass is mpp < 160 meV, which arises from the recently quoted EXO-200 experiment |62j 
T°J'^^(^^®Xe) > 1.6 X 10^® yr at 90% C.L. This limit will be improved within a not too distant future. The GERDA 
“phase-H”experiment [63l[64] is expected to reach iJ^Ge) > 2 x 10^® yr, corresponding to W/j/g < 100 meV. A 
bolometric CHORE experiment, using |65| . is currently under construction. It has an estimated sensitivity 

around ~ 10^® yr, which corresponds to mpp < 50 meV. There are also proposals for ton-scale next-to- 




next generation Oi//3/3 experiments with ^®®Xe [551123 and ^®Ge [531155] claiming sensitivities over ^ 10^^ yr, 

corresponding to ^12 — 30 meV. For recent reviews, see for example Ref. [69j . Consequently, as it can be seen 
from Eq. (36), our model predicts Ty 2 ^ at the level of sensitivities of the next generation or next-to-next generation 
Oiypp experiments. 


IV. QUARK MASSES AND MIXING 


Using Eq. (15), together with the product rules for the Ty group given in Appendix]^ and considering that the 
components of the SU{5) singlet Tj scalar triplet ^ acquire the same VEV as shown by Eq. (10), which is set to be 
equal to AA as the VEV of the Z 12 charged scalar cr as indicated by Eq. © (A being the cutoff of our model), we 
find the quark mass matrices: 


Mrr = 


Oil A 

(&),3 

13 ^ 


a 


(&) a 4 3).2 

22 ^ ^23 ^ 


a. 


a. 




23 


33 




(37) 




0 

0 


0 


0 

0 


3 ') = diag{md,rns,rnb) {v^j)j 


(38) 
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where A = 0.225 and the 0(1) dimensionless couplings in Eqs. (371 and (38) are given by: 


*12 


= - 


(4) 

= 2^2(712+721) 

( 2 ) 
v\ ^ 


= 2 ^ 2(723 + 732 ) 

1 


V 




2/3iU$^ 


*11 




,( 2 ) 


= 4^2711^, 


V 

,(3) 


*13 


.,(3) 


= 2 V^( 7 i 3 + 731) 


= 4 V 2722 —, 


„(4) 


(D) 4 

a\ — - 


{a2V%^ + 2/32 


43'=4v^733^, 


.(4?) 


CLo — 


2/33U$) . (39) 


From Eq. (38) it follows that the CKM quark mixing mixing matrix does not receive contributions from the down-type 


quark sector, meaning that quark mixing arises solely from the up-type quark sector. Consequently, having a realistic 
up-type quark masses and quark mixing angles requires that the mass matrix for up-type quarks to be given by: 


Mu = VcKMdiag (to„, mc,mt) V) 


CKM- 


(40) 


Using the standard parametrization for the CKM quark mixing matrix, together with the relations rriu = oA^mt and 
rxxc = bX^mt, where a and b are 0(1) coefficients, we get the mass matrix for up-type quarks described in Eq. (37), 
with entries exhibiting different scalings in terms of powers of the Wolfenstein parameter A = 0.225. Thus, from the 
requirement of realistic up quark masses and quark mixing angles with 0(1) dimensionless couplings in Eq. (37), we 


hnd a A suppression in the 11 entry of the up-quark mass matrix instead of a A one. We have numerically checked 
that a A® supression in the 11 entry of the up-quark mass matrix, with 0 ( 1 ) dimensionless couplings, is consistent 
with realistic up-quark masses and quark mixing angles. 

Assuming that the quark mass and mixing pattern is caused by the breaking of the Z 3 , Z 4 and Z 12 symmetries, 
to simplify our analysis, we adopt a benchmark where the dimensionless charged fermion Yukawa couplings are 
approximately equal. Specifically, we set 


7ii = (1-y 


1/2 


7ie*'^i, 


7i 3 = 731 = 72 ( 1 - y 


v2 \ 3 




7i 2 = 721 = -7ie*‘^U 


723 ~ 732 — 72) 


( 1 ^ 

722 = 7i 1 1 - y 


-1/2 


* = 1,2,3, 


(41) 


with 7 ]^, 72 , ai and /3j (/ = 1,2,3) real 0(1) parameters. Our benchmark of nearly equal charged fermion Yukawa 
couplings given in Eq. (41), which we have numerically checked is consistent with realistic up-quark masses and quark 


mixing angles, is also adopted in order to reduce the number of free effective parameters in the quark sector of our 
model. There is no tuning in the parameters of our model. Let us note that the exactly equal dimensionless quark 
Yukawa couplings do not allow generating the up and charm quark masses. In Appendix]^ we give another possible 
motivation for the approximate universality of dimensionless couplings, which can be studied beyond the present 
model, adding new symmetries. 

Besides that, for simplicity we assume that the complex phase responsible for CP violation in the quark sector arises 
solely from the up-type quark sector, as indicated by Eq. (411. In addition, to simplify the analysis, we fix = 1, 
as suggested by naturalness arguments. Consequently, the up-type quark mass matrix reads 




Mu = 


('W) 

AU)- 
*1 


1/2 




-.(u) 


X^e^4>2 


AU) 


V 


AU) 


(‘-^) 


5„i02 
3 


AU) 


A^e*‘^3 


('W) 


- 1/2 


('W) 


l-4rl A^e*'^3 


\ 


2 

(U) 

2 

1 


V 

71’ 


(42) 




As seen from the above formulas, the quark sector of our model contains ten parameters, i.e. A, a'l^\ **7^, 

**7^ ^ 4 id the phases </); (^ = 1,2, 3), to describe the quark mass and mixing pattern, which is characterized by 
ten physical observables, i.e., the six quark masses, the three mixing angles and the CP phase. Out of the ten model 
parameters two of them, A and are fixed, whereas the remaining eight are htted to reproduce the six quark 

masses and four quark mixing parameters. In Table m we show the experimental values of the physical observables 
in the quark sector, together with our results obtained for the following best-fit values of the model parameters: 


*7^ 

AD) 


~ 1.96, 


AU) 


~ 0.74, 


ai ' ~ 2.54, 


^ 0.58, 


61 ~ 10.94°, 
~ 1.42. 


62 6.02°, 


~ 21.65°, 


(43) 
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Observable 

Model value 

Experimental value 

mn{MeV) 

0.86 

-L-^0 o.45 

mc{MeV) 

673 

635 ± 86 

mt{GeV) 

174.2 

172.1 ±0.6 ±0.9 

md{MeV) 

2.9 

2 

^•^-0.4 

nisiMeV) 

57.7 

^'•'-15.7 

miiGeV) 

2.82 

c\ qo+0.09 

Z.OZ Q Q4 

\Vud\ 

0.974 

0.97427 ±0.00015 

|U„e| 

0.2257 

0.22534 ± 0.00065 

iKi, 

0.00305 

0 nnQc:i +0-00015 

U.UUO01 Q QQQ;^4 


0.2256 

0.22520 ±0.00065 


0.97347 

0.97344 ± 0.00016 

\Vcb\ 

0.0384 

0.04121Hm5 

\vj 

0.00785 

0.00867l°;ro3? 

Usl 

0.0377 

U.U^U^0.0005 

\Vtb\ 

0.999145 

'^•^^^-‘-^^-0.000046 

J 

2.32 X 10"® 

X 10"® 

5 

64° 

68 ° 


Table III: Model and experimental values of the quark masses and CKM parameters. 


We use the experimental values for the quark masses at the Mz scale, reported in Ref. ([70]) (which are similar to 
those in d]), whereas the experimental values of the CKM matrix elements, the Jarlskog invariant J and the CP- 
violating phase 5 are taken from Ref. [5] . Let us note that the agreement of our model with the experimental data is 
as good as in the models of Refs. [72U79] . and better than in Refs. [8T1I80H87] . The following comparison of our model 
with these models could be in order. Despite the similar quality of the data description our model is more predictive 
than the model of Ref. since the latter, focused only on the quark sector, has a total of 12 free parameters, 

whereas the quark sector of our model is described by 8 free effective parameters that are adjusted to reproduce the 
10 physical observables of the quark sector. The models of Ref. [75], Ref. [7^ l80] . Ref. [82], Refs.[79l HT] |87|, Refs. 
[TlIMZilHS], Refs. (Huns] and Ref. [53] possess in the quark sector 6, 7, 8, 9, 10, 12, and 13 free parameters. The 
total number of the effective free parameters of our Ty flavor adjoint SU (5) GUT model is 16, from which 2 are fixed 
and 14 are fitted to reproduce the experimental values of 18 observables in the quark and lepton sectors. On the other 
hand the SU{5) model with T' 0 Zyi ® symmetry of Ref. [3T] has nine parameters in the Yukawa sector for the 
charged fermions and the neutrinos. However it does not account for CP violation in the quark sector, whereas our 
model does. Furthermore, the values of the physical observables we derived in our model for both quark and lepton 
sectors exhibit a significantly better agreement with their corresponding experimental values than those derived in 
Ref. [3T] within the SU (5) model with T' ® Zyi ® Z'^^ symmetry. 


V. GAUGE COUPLING 5U(5) UNIFICATION 


In the previous sections, we analyzed the possibility of describing the quark and lepton masses and flavor mixing 
within the framework of our model. This analysis was based on the symmetries of the model and particular field 
assignments to the symmetry group representations, which allowed us to single out several effective parameters 
completely determining the lepton and quark mass matrices. It is a notable property of the model that the SM 
fermion mass matrices depend on the “fundamental” parameters of the model Lagrangian only through this set of a 
few effective parameters. Now we turn to more subtle aspects of the model depending on the details of the non-SM 
components of the SU{5) multiplets, as well as on the “fundamental” parameters, which may have crucial impact on 
its ultraviolet behavior. 

As is well known, there are many extensions of the SM which lead to gauge coupling unification (GCU) and also 
successfully fulfill all the constraints coming from fermion masses, proton decay, and perturbativity. In particular, 
models based on supersymmetric and also nonsupersymmetric (non-SUSY) SU{5) unification have been widely studied 
in the literature (22 [55]. For non-SUSY SU{5) scenarios, the unification of gauge couplings can be as good as or 
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better than in the MSSM, despite the fact that the number of extra fields beyond the SM is smaller. More restrictive 
conditions such as the possibility of implementation of an appropriate neutrino mass generation mechanism, and 
compatibility with the existing phenomenological, cosmological and astrophysical constraints, require some specific 
properties for this extra field content. As was already pointed out in Sec. II, the SU{5) scalar representations with 
the minimal number of Higgs bosons needed to generate the fermions masses and mixings are one 24s, one 45s, four 
5s ’s, and twelve Is’s. If an extra fermionic 24f representation is also included, a simple configuration of the extra 
fields allowing the type-I and type-III seesaw mechanisms of neutrino mass generation can be constructed with masses 
below the GUT scale. As will be shown later the condition of the GCU and compatibility with the lower experimental 
limit on the proton decay half-life predicts the masses for the extra scalar fields within the LHC reach. This opens 
up the possibility for experimental tests of the considered models. On the other hand, the dark matter constraints on 
the scalar sector will lead us to the conclusion that the masses for the type-I and type-III fermionic seesaw mediators, 
mp^ and nip^, should be at high energies far beyond the TeV ballpark. 

It is noteworthy to mention that our configurations rely on fine-tuning which separate light from heavy degrees of 
freedom. However, this issue is an unavoidable problem of grand unified theories. Evidence of this is the standard 
SU{5), which suffers from what is known as the doublet-triplet splitting problem i.e, the fact that the SM Higgs 
doublet has a mass of rrih ^ 125 GeV while the colored triplet in the 5 irrep must have a mass of the order of the 
GUT scale in order to prevent proton decay. Although several solutions to this problem have been suggested in the 
literature (for a short review, see Ref. [89j 1. we do not consider here any particular one. Instead, we view this as a 
fine-tuning problem and, in fact, to make the exotic particles in our model light will require, in general, additional 
fine-tunings. We assured ourselves that the large number of free and uncorrelated parameters in the scalar potential 
will allow us to reproduce the required mass differences. There are non-SUSY SU(5) models which also heavily rely 
on fine-tuning among states, and which do not consider formal solutions to this problem, as seen in Refs. [301 sg. 
We are aware of this fine- tuning problem, whose formal solution goes beyond the scope of this work and is deferred 
for a future publication. 


1. Setup 


The scalar sector of our model is composed, as described in the previous sections, by the 5s, 24^, and 45^ irreps of 
SU{5) with the following SU{3)c x SU{2)l x U{ 1 )y decompositions: 

5s = (1, 2 )i/2 © (3, 

24s = (I, l)o © (Ij 3)o © (8, l)o © (3, 2 )_ 5/6 © (3, 2 ) 5 / 0 , 

= .nil © .^3 © .rig © “(3.2) © ^( 3 , 2 )j 

45s = (1, 2 ) 1/2 © (3, l)-i/3 © (3, 3 )_i/3 © (3, 1 ) 4/3 © (3,2)_7/6 © (6,1)-i/3 © (8, 2 ) 1 / 2 , 

= $1 © $2 © 4*3 © $4 © $5 © ^6 © $7- 

(44) 


The SM Higgs doublet is embedded in the 5s. The adjoint 24s representation triggers the breaking of SU{5) to the SM 
at the GUT scale. As we previously commented, this particle content should be further extended in order to account 
for the nonzero neutrino masses. We introduced a fermionic 24/ irrep of the SU{5). It has the SM decomposition 

1)0 © (li 3)0 © ( 8 ,1)0 © ( 3 , 2 )_ 5/0 © ( 3 , 2)5/0, 

= Po ® P 3 ® Ps ® P(3,2) ® P(3,2)- (45) 

In this case the SU (5) singlet Majorana neutrino with the heavy SM singlet Pq mediate type-I seesaw mechanism 
while the SU{2)l triplet gives rise to the type-III seesaw mechanism. Combining these three extra fermions, Nji, 
Pq, and p^, with some of the scalar fields from the 5s, 24s, and 45s irreps, we construct the simplest benchmark 
configurations that unify the gauge couplings within the SU{5) and satisfy some general requirements in order to 
guarantee their phenomenological viability. These requirements are 

(i) Perturbative SU{5) unification: This means that gauge couplings unify as well as or even better than in the MSSM, 
and the value of ac is in the perturbative regime. Note that we are not necessarily imposing the exact unification of 
the gauge couplings at the GUT scale (me)- Rather, we allow for a difference of the gauge couplings at mo falling 
into the area of the MSSM ’’nonunification triangle” [soiini]. 

(ii) Proton decay: There are some specific fields contributing to proton decay. The dimension-six proton decay 
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operators are mediated by the superheavy gauge bosons, usually named leptoquarks, in the 24 irrep: P( 3 ^ 2 ) ® P (3 2 ) = 
(3, 2, —5/6)0(3,2, 5/6), which must be heavier than 3x 10^® GeV to satisfy the experimental lower bound on the proton 
decay lifetime. Here, we assume that these fields live at the GUT scale. In addition, we deal with field configurations 
which, in almost all the parameter space, fulfill the current constraint from Tp_>^ 0 e+ > 10^^ years [92l[93]. This, 
through the relation for the proton decay half-life, T = leads to a GUT scale of me ^ 3 x 10^® GeV. 

(iii) Fermion masses: In particular, neutrino mass generation through the type-I and type-III seesaw mechanisms 
[miMj. The configurations should then contain at least one copy of the fermionic fields Nn, pg, and p^, as described 
before. 

(iv) Nontrivial phenomenology: Among the models passed through the above three conditions, we select those models 
which may provide a distinguishing signal at the LHG. 

In what follows, we analyze some of the ’’minimal” benchmark models, which lead to correct unification and also 
fulfill the above listed requirements, with as simple a field content as possible. As the preferable models, we consider 
those which are phenomenologically rich enough, in the sense that some of the fields (being colored) could give rise 
to certain resonances at the LHC. The analyzed models will lead to an available parameter space for the different 
masses of the scalars and the type-I, type-III seesaw fermionic mediators. 


2. The Models 


For simplicity, and following the notations of Ref. [33] , we rewrite the masses of the fermionic 24 j given in 
as follows: 


Eq. (19) 


mpj = m — 3e*“A , 
mpg = m -I- 2e*“A , 

”^^3.2) =""'>(3.2) =w-l/2e“A', 

rupg = m — e*“A . (46) 


where a is the relative phase between yi and y^. For the particular case where a = 0, the parameters A and m are 
then defined as m = y 2 Vri and A = yavs.^ (where ys = ys/V^)- 

In order to deal with the simplest models, we look for configurations where the only contribution from the fermionic 
24/ to the RGE flow comes from the fermionic type-I and type-III seesaw mediators. We assume the other 24/ 
components, pg, P( 3 p), and P (3 2 )j bave no RGE effect, being as heavy as the GUT scale. It is worth reiterating that 
in the analyzed benchmark models the GCU is achieved having a few particles with masses below the GUT scale: the 
fermions Nji, pg, pg plus some of the extra scalar fields from the 5^, 24^, and 45^ multiplets. This kind of spectra 
can be easily obtained by fine-tuning Eq. (46), imposing that the mass of the remaining fermionic fields in the 24/ 
lives at the GUT scale or above. 


Searching for the models, we keep our analysis at the one-loop level. It could be easily extended to two loops. However, 
this sophistication makes no impact on our final results and conclusions. 

The master equation for the running of the inverse gauge couplings at the one-loop level is 


a* 



27 T 



with z = 1, 2, 3 


(47) 


The effective one-loop RGE coefficients, taking into account the thresholds from particles with masses m/, are given 
by 


= ^0{T-mf)b{. (48) 

/ 

The contribution of each particle b{ is calculated according to 

h = -yr,(RG) + \%{Rf) + \r{Rb). (49) 

where T{Rj) are the Dynkin indexes of the representations Rj to which belong I = G, F, and H-the gauge bosons, 
fermions, and scalars, respectively. They are defined as T{R)Smn = Tr(Tm(i?)T„(i?)), with Tm{R) being generators 
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in the representation R. For the lowest-dimension representations of SU{N) they are T{fundamental) = 1/2, 
T{Adj) = N. For the SM we have {b§^,b 2 ^,bf^) = (-7,-19/6,41/10), which correspond to the contributions 
of the SM fermions and one copy of the SU{2)l Higgs doublet. The additional non-SM fields introduce extra 
contributions Ab^ to these coefficients bi = bf^^ + Abj. In Table IV Appendix the Ab^ contributions of the fields 


in the 5s, 10s, 24s, 24/, and 45s SU{5) repres entat ions are shown. In order to correctly unify the gauge couplings 
and fulfill all the requirements (i)-(iv) in Sec. V 1 which we impose on the models, these Ab^ coefficients should 


obey certain conditions. Some of the simplest configurations of the non-SM fields with the masses below the GUT 
scale, which obey these conditions, are listed in Table [V] AppendixAll the listed field configurations have a highly 
split mass spectrum with type-I and type-III seesaw mediators at ms ~ 10^^ GeV and remaining scalars at m^p = 2 
TeV. Let us note that since among the scalars there are the color octets, we use 2 TeV as a limit recently established 
on the mass of color octets by the GMS and ATLAS collaborations from the dijet pair signature searches. We 
have purposely chosen the latter scale to be low enough so that the colored scalars are within the LHC’s mass reach, 
while the large value for the seesaw scale ms is in agreement with the small neutrino masses. Let us consider two 
simplest models (1) and (2) from Table [v| 


Model (1) This is the simplest of all the benchmark models passing our conditions (i)-(iv) in Sec. V 1[ including the 
unification of the gauge couplings and smallness of neutrino masses. However, if we fix the masses of the scalar fields 
at TTiArp = 2 TeV and the fermionic seesaw mediators pg 3, Np at ms = 10^^ GeV, as it is done for all the models 
in Table [v| we get the GUT scale me = 2.7 x lO^^GeV, which is in slight tension with the limit imposed by the proton 
stability toq > 3 x 10^^ GeV. This flaw can easily be cured by allowing the masses of the fields to vary independently 
in the range 2 TeV< m^ijmss ^ ws < mo- Now, as shown in Fig. a significant part of the model parameter 
space corresponds to the GUT scale in the range allowed by the proton decay constraint. From the left panel of this 
figure we can see that in this part of the parameter space the masses ms of the seesaw mediators pg 3, Np could be in 
the ballpark of [10®•^ GeV. The right panel of Fig. 0 shows that in the part of the parameter space consistent 

with the proton decay constraints, the mass of the scalar octet Ss is relatively low, 2 TeV< msg ^ 10® GeV. Thus, 
there is a chance for Sg to be within the mass reach of the current run of the LHG. However, its production cross 
section suffers from several suppression factors. The process of single-Sg production is only possible in the gluon 
fusion gg —>■ Sg via a loop with two or three internal Sg, which is suppressed by the large m^g. On the other hand, 
as seen from Eq. (15), it cannot be produced in qq ^ 2g. The tree-level pair production process gg —>■ 5gSg, being 
not suppressed in the amplitude, has a high threshold of 2mHg. If produced, 5g decays in a unique channel Sg —>■ gg. 
However, the corresponding signal at the LHG could be challenging from the viewpoint of the identification of its 
origin. Below we consider another benchmark model with a more distinctive signature of the color octet from 45. 


A viable dark matter (DM) candidate in non-SUSY models is an issue requiring special efforts. References [43l [96] 
recently proposed the neutral component S® of the S 3 ^ (1, 3 )g C 24^ as a scalar cold dark matter (GDM) candidate. 
The necessary condition for its stability is the vanishing of the trilinear coupling of with the SM Higgs dou¬ 

blet, H. Also, its VEV must be zero, (S®) = 0. These conditions, as shown in Ref. [43], can be implemented by an 
ad hoc fine-tuning of the scalar potential parameters. Unfortunately, no proper custodial symmetry, protecting this 
fine-tuning, can be incorporated in the adjoint SU{5) framework [43] . Despite this complication, we study constraints 
on our model from the assumption that the S® is a GDM candidate. It has been shown in Ref. [35] that the thermal 
relic abundance could be compatible with the observed DM abundance = 0.110 ± 0.005, if ~ 2.5 TeV. 

We use this fixed value as a condition allowing the presence of a viable GDM particle candidate S® in the model. In 
Fig-i we show limitations on the model parameter space, taking into account this condition. As seen, the lowest 
bound on the seesaw mediator mass is ms > 10^® GeV. This large value of the seesaw mass scale perfectly accounts 
for the smallness of the neutrino masses. 


Model (2): This is the next simplest benchmark model from Table It contains a scalar color octet 
4)7 ^ ( 8 , 2 ) 1/2 C 45. Its production and possible signals at the LHC have been studied in Ref. [97] (for the ear¬ 
lier studies of the color octet scalars at the LHC see, for instance. Ref. [M])- The neutral and charged components 
Sp j, of this multiplet, unlike the Sg in Model (1), can be singly produced at tree level in quark-antiquark annihi¬ 
lation qq —)■ j, q'q —> and in the gluon fusion at one-loop level via &-quark loop, which is much less suppressed 
than the loop-induced production of a single Sg. The tree-level pair production, dominated by the gluon fusion 
gg —)■ S^S^yS^S^, is approximately an order of magnitude smaller than the single-S” production at ms ~ 2 TeV, 
according to Ref. IS?]. However, the single-S* production cannot overcome the SM background with all the versatile 
kinematical cuts applied in m, and therefore has no observational prospects. On the other hand, it has been shown 
that the S'-pair production, even having a smaller production cross section, may get larger than the background in 
the 46-tagged final-state jets with a cut Pp > 800 GeV. 

Repeating the analysis made for the previous model, we find masses m^j compatible with the proton decay constraints 
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me > 3 X lO^^GeV. As shown in Fig. we have 2 TeV< m<^,^ < 10^^ GeV. 


Figure 1: The parameter space passing the conditions (i)-(iv) in Sec. V 1 The region compatible with the proton decay 


constraints ma > 3 x 10'^°GeV is explicitly shown in light blue. The upper two panels correspond to Model (1) and the lower 
one to Model (2) considered in the text. In the latter case, the mass of the extra scalar electroweak doublet >l>i is fixed at 2 
TeV. 




2 4 6 8 10 12 14 


Login (m<I>7/GeV) 


Let US note also that, even if we fix the mass of all the scalars at 2TeV and the seesaw mediators at GeV, as is 
done for the models in Table the GUT scale is still high enough to be safe, facing new possible improvement of 
the proton decay constraints in the future experiments. Now, as with Model (1), we consider the constraints imposed 
on the present model by the interpretation of the S 3 as a CDM candidate. As we discussed above, this requires that 
2.5 TeV. Scanning the parameter space, as with Model (1), we find the plot shown in the right panel of Fig. 
from which we derive the lower bound ms ^ 10^^ GeV. 



LogjoCms/GeV) 


Logio(ms/GeV) 


Figure 2: Dark matter constraints: allowed parameter space in ms and mss- The left and right panels correspond to Models 
(1) and (2), respectively. 
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VI. CONCLUSIONS AND DISCUSSION 


We proposed a version of the adjoint SU{5) grand unification model with an extra T 7 0 Z 2 0 Z 3 0 Z 4 0 Z 4 0 Z \2 
flavor symmetry, which successfully describes the SM fermion mass and mixing pattern. The model has in total 
16 effective free parameters, from which 2 are fixed and 14 are fitted to reproduce the experimental values of 18 
observables in the quark and lepton sectors i.e., 9 charged fermion masses, 2 neutrino mass-squared splittings, 3 
lepton mixing parameters, 3 quark mixing angles, and 1 CP-violating phase of the CKM quark mixing matrix. The 
observed hierarchy of charged fermion masses and quark mixing angles is a consequence of the 03 004 0 Zi 2 symmetry 
breaking, triggered by the SU{5) scalar singlets cr, r, and tp, charged under this symmetry, and which acquire VEVs 
at very high-energy scale, close to the GUT scale. The non-SM fermion spectrum of our model is composed of a single 
heavy SU (5) singlet right-handed Majorana neutrino Nji and two more fermionic fields: an electroweak singlet pg 
and triplet pg, both from the adjoint 24 irrep of SU{5). Thus the light neutrino masses arise in our model from type-I 
and type-III seesaw mechanisms mediated by these fields. The smallness of neutrino masses is a consequence of their 
inverse scaling with respect to the masses of these three seesaw mediators as well as the quadratic proportionality to 
the presumably small neutrino Yukawa couplings. The model predictions for the physical parameters in the quark 
and lepton sectors are in good agreement with the experimental data. The experimentally observed deviation from 
the trimaximal pattern is implemented by introducing two T-j triplet scalars x s-iid singlets under SU{5). The 
model predicts an effective Majorana neutrino mass, relevant for neutrinoless double beta decay, with values mpp = 
4 and 50 meV, for the normal and the inverted neutrino spectrum, respectively. In the latter case, our prediction 
is within the declared reach of the next-generation bolometric CUORE experiment [5S] or, more realistically, of the 
next-to-next-generation tone-scale 0://3/3-decay experiments. 

According to our model, the leptonic Dirac CP-violating phase is vanishing. In view of the current experimental 
trend, it could get into tension with the observations. The latest fit of the neutrino oscillation experimental data by 
the Valencia group. Ref. [5], indicates a nonvanishing CP-violating phase at the Itr (less than 2cr) level for the normal 
(inverted) neutrino mass hierarchy. The Bali group in their fits claims slightly larger significance, but also below the 
3 ct level [5S]. Thus, with the current experimental significance, essentially below than the golden 5cr, our model is 
not yet in trouble with the CP. On the other hand, we are aware of the necessity of reconsidering this aspect of our 
approach, which will be done elsewhere. 

In the last section of this paper, we studied the compatibility of our model with certain physical conditions expected 
from a plausible GUT model. Among them we considered the gauge coupling unification, the proton stability con¬ 
straints, the smallness of the active neutrino masses, and nontrivial LHC phenomenology. Towards this end, we 
specified the simplest benchmark models of the non-SM particle spectrum lighter than the GUT scale and which meet 
these conditions. We examined two of them and found that they may give rise to observable signals at the LHC, as 
well as contain a viable scalar CDM particle candidate. In order to account for the observed DM relic abundance in 
the Universe, the latter imposes certain constraints on the model parameter space. From this DM condition we found, 
in particular, the lower limits in the mass of the seesaw type-III mediator (ms), for the two simple benchmark models 
analyzed in the paper. It is worth noting that Model (1) is manifestly falsifiable, as seen from Figj^ A relatively 
small improvement of the proton decay lifetime lower limit up to ^ 6 x 10^® GeV would reject the model. 

Various aspects of the adjoint SU{5) scenario have been studied in the literature [HI [43l llOOj . In Refs. [HI [43], 
the fact that is the lightest field in the 24 offers a viable way to understand the baryogenesis via leptogenesis. In 
this work, we also considered the field as the lightest component of the 24, and therefore, the baryogenesis can be 
accomplished in our model in the same way. Note also that in this mechanism uno], the decays of the p^ create a 
lepton asymmetry, which then is converted in a baryon asymmetry by the sphalerons. Imposing the condition of the 
successful leptogenesis and for the normal hierarchy of neutrinos, one finds that the mass of should be large, which 
is in agreement with the large p^ mass values predicted in the present paper. 

Finally, as was previously noticed in Ref. [13], the mass of pg must be heavier than 10® GeV—10^ GeV in order to 
satisfy the constraints from the big bang nucleosynthesis for the GUT scales larger than 3 x 10^® GeV. This condition 
is also consistent with our results, where the mass of this field is set around the GUT scale. 
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Appendix A: The product rules for T7 

The group Ty is the minimal non-Abelian discrete group having a complex triplet. The discrete group T 7 has 21 
elements and 5 irreps, i.e., one triplet 3, one antitriplet 3, and three singlets Iq, li and I 2 |12| . Furthermore, the 
T 7 group is a subgroup of SU{3) and A{3N^) with N = 7 and is isomorphic to Z 7 x Z 3 . The triplet and antitriplet 
irreducible representations can be defined as follows [T^ : 


3 = 



3 = 



(Al) 


The triplet and antitriplet T 7 tensor irreducible representations satisfy the following product rules: 




I ^ 0:22/5+w^''a; 4 y 3 )ifc- 

fc=0,l,2 


(A2) 


(A3) 


(A4) 


Whereas the tensor products between singlets are given by the relations 

(a:)io( 2 /)io = {x)u{y)i2 = = {xy)-i 

{x)iAy)ii = {xy)i2, 

{x)i2{y)i2 = {xy)i,, 

the product rules between triplets and singlets read 


(A5) 



yxi(e) 
yX2{5) 

/ 3(3) 


(A 6 ) 


where uj = . The representation Iq is trivial, while the nontrivial li and I 2 are complex conjugate to each other. 

Some reviews of discrete symmetries in particle physics can be found in Refs. pTlIT^ISQlllOlj . 


Appendix B: Ou the universality of Yukawa couplings 


The aforementioned scheme of approximate universality of dimensionless couplings can be justified by adding an extra 
^24 symmetry and four SU{5) scalar singlets, assigned as T 7 trivial singlets, as well as by setting their VEVs to be 
equal to AA, with A = 0.225, one of the Wolfenstein parameters, and A being our model cutoff. One of the four SU{5) 
scalar singlets, namely S, can be assumed to have the same Z 12 charge as the scalar field a of our model and can also 
be made charged under the new Z 24 symmetry and neutral under the remaining cyclic symmetries. The remaining 
three SU{5) scalar singlets, namely {i = 1,2, 3), can be assumed to be only charged under the new Z 24 symmetry. 
These four new scalar fields will transform under the Z 24 symmetry as follows: 

S->■ S, Ai^-fAi, A2->e*S'A2, A 3 ->■ A 3 . 


(Bl) 
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The aforementioned Z 24 charge assignments will generate the following Z 24 neutral combinations of the new scalar 
fields: 


S^A^S^Al 


S^Al 


(B2) 


These Z 24 neutral combinations will give rise to the following Yukawa operators invariant under the symmetries of 
the model: 

ijfcip,T,(l)^(2),T,(l) ^1 zifcip,T,(2)^(3),T,(2) ^2 ijfcip,T,(l)^(3),T,(3) ^3 z4fcZp,T,(3) ^(3) ^3 /rio'i 

ij P kl ^8 ’ b P kl ^6 ’ ij P kl ^5 ’ ^ ^ ij ^^p ^ kl ^5 ' 

The first two aforementioned Yukawa operators will contribute to the 11 and 22 entries of the up-type quark mass 
matrix, whereas the last two operators will contribute to the 13 and 31 entries. These new contributions will be 
proportional to A®, A®, and A®, respectively, and thus will introduce deviations from the exact universality in the 
dimensionless Yukawa couplings. These aforementioned contributions will correspond to the first-order term in the A 
expansion of the expressions given in Eq. (41). 


Appendix C: Simple benchmark models 


IV 


In 


The contributions to the Ab coefhcients for each held in the 5s, 24s, 45s, and 24/ reps are shown in Table 
Table |Vj the conhgurations which generate neutrino mass in agreement with the mechanism described in Sec. Ill 
are shown. Each particle content, added to the SM, lead to ”SM+X” conhgurations unifying gauge couplings almost 
equal to or better than the MSSM, i.e, each one of the simplest conhgurations satishes: ^ 


<^ 2 MSSMi'^G) - o^imssm(.'^g) and 3 X 10^® GeV < toq < 
the actual experimental bounds. 


10^® GeV in order to obtain proton lifetimes allowed by 


SU{5} Rep 

Fields 

^^1,2,3 

5s 

Hi,H2 

(i,l,0),(i,0.1) 

24s 

Hi, H 3 , Hs, S(3,2), ^(3,2) 

(0,0,0), (0,1,0), (0,0,1), (1,1,1), (1.1,1) 

45s 

^1) ^2? ^3? ^45 ^55 *^6) ^7 

(T.^o). (T.0.1) . (1,2,1) (A,0.1), (A.1.1) (A,0.1) (4.|.2) 

24/ 

Pi? Ps’ PS’ P(3,2)) P(3,2) 

(0,0,0), (0,1,0), (0,0,2), (1.1,|), (1,1,1) 


Table IV: Contribution to the running of the gauge couplings of the fields from the 5s, 24s, 45s and 24/ representations of 
the SU{5). The seesaw particles mediators Nji and Pq 3 are added at the 10^^ GeV, contributing to the running of the gauge 
couplings with (0, 4/3, 0) in the range [10^‘^GeV, ma]- The contribution of the fields in the 24/ is twice the contribution of 24s 
owing its fermionic nature. 
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Table V: Simple X configurations, leading to the GCU and generate the neutrino mass via type-I and type-III seesaw mech¬ 
anisms. The masses of all the extra scalar particles and Hi are about mjvp = 2 x lO^GeV while the the fermionic seesaw 
mediators Nii,pQ 3 , are much heavier with the masses at the scale ms = 10^"^ GeV. The A 6 s RGE coefficients correspond to 
the contribution of the scalars. The seesaw mediators contribute with an extra (0,4/3,0) above the scale lO^^GeV. 
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